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Gravitational potential of a homogeneous circular torus: 
new approach 
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ABSTRACT 

The integral expression for gravitational potential of a homogeneous circular torus 
composed of infinitely thin rings is obtained. Approximate expressions for torus po- 
tential in the outer and inner regions are found. In the outer region a torus potential 
is shown to be approximately equal to that of an infinitely thin ring of the same mass; 
it is valid up to the surface of the torus. It is shown in a first approximation, that 
the inner potential of the torus (inside a torus body) is a quadratic function of coor- 
dinates. The method of sewing together the inner and outer potentials is proposed. 
This method provided a continuous approximate solution for the potential and its 
derivatives, working throughout the region. 
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1 INTRODUCTION 

Toroidal structures are now detected in astrophysical objects 
of various types. Such objects are, for example, the ring 
galaxies, where a ring of stars is observed. In some galaxies, 
the ring-like distribution of stars is believed to be due to 
collisions of galaxies, as, for example, in M31 (Block et al. 
1987), and Arp 147 (Gerber et al. 1992). The analysis of 
the SDSS data (Ibata et al. 2003) indicates the existence 
of a star ring in the Milky Way on scales of about 15-20 
kpc, which is believed to be originated from the capture 
of a dwarf galaxy. Obscuring tori are observed in central 
regions of active galactic nuclei (AGN) (Jaffe et al. 2004) 
and play an essential role in the unified scheme (Antonucci 
1993; Urry & Padovani 1995). Ring-like structures exhibit 
themselves in dark matter as well. An example can be the 
galaxy cluster C10024+17 where a ring-like structure has 
been found in distribution of dark matter with the use of 
gravitational lensing method (Jee et al. 2007). In the Milky 
Way, the rotation curves together with the EGRET data can 
be explained by existence of two rings of dark matter located 
at distances of about 4 kpc and 14 kpc from the Galaxy 
center (de Boer et al. 2005). Such toroidal structures can 
possess a significant mass, and thus gravitationally affect 
the matter motion. 



tational potential 0. Kondratyev (2003) has returned to this 
problem for the first time. In this work an exact expres- 
sion for the potential of a homogeneous torus on the axis of 
symmetry was obtained. In (Kondratyev 2007) the integral 
expressions for a homogeneous torus potential were found 
using a disk as a primordial gravitating element. Stacking up 
such disks will result in a torus with the potential equal to a 
sum of potentials of component disks. However, it is evident 
that any integral expressions are problematic to use both 
in analytic studies and in numerical integration of motion 
equations, and also in solving the problems of gravitational 
lensing. B. Kondratyev et al. (2009, 2010) have obtained an 
expansion of torus potential in Laplace series, but showed, 
however, that such an expansion is impossible inside some 
spherical shell. 

In this paper we propose a new approach to investiga- 
tion of the gravitational potential of a torus. Special atten- 
tion has been paid to finding approximate expressions for the 
potential, which would simplify investigation of astrophys- 
ical objects with gravitating tori as structural elements. In 
contrast to (Kondratyev 2007), we used an infinitely thin 
ring as a torus component. Such ring is actually a realiza- 
tion of a torus, with its minor radius tending to zero and 
the major one equaling the ring radius. Using such an ap- 
proach, we obtained an integral expression for the potential 



B. Riemann devoted one of his last works to the grav- 
itational potential of a homogeneous torus (see Collected 
papers, 1948). This work remained unfinished however. For 
over a century, no attention has been paid to a torus gravi- 



1 In electrostatics, a potential of conducting torus shell is used 
(Smythe 1950) that is much easier than the case when torus den- 
sity is uniform inside its volume. 
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Figure 1. On the top: a 3D scheme of a torus configuration; on 
the bottom: schematic cross-section of a torus as a set of infinitely 
thin rings. 



of a homogeneous circular torus (Section 2) and approxi- 
mate expressions for the potential in the outer (Section 3) 
and inner (Section 4) regions. In Section 5, the method of 
determining a torus potential for the entire region is sug- 
gested. 



2 GRAVITATIONAL POTENTIAL OF A 
HOMOGENEOUS TORUS 

Compose a torus with mass M, outer (major) radius R and 
minor radius Rq , of a set of infinitely thin rings - component 
rings hereafter, - (see Fig. 1), with their planes being parallel 
to the torus symmetry plane. Select a central ring with mass 
M c and radius 7? from a set of rings composing a torus. The 
potential produced by this ring at an arbitrary point P(x, z) 
has a form: 

„ „ , , GM C , / x z \ 
where the dimensionless potential of the infinitely thin ring 



(2) 



— is the complete elliptical integral 

nete: 
AxR 



-m sin 2 8 

of the first kind, and its parameter is 



(x + R) 2 + z 2 



(3) 



The potential at a point P(x, z) produced by an arbitrary 
ring with radius R' and mass M r located in the torus at a 



GM 

<p r (x,z;M r ,z ) = 



(4) 



where the expression for (f> r is obtained by substitution in 
© of a kind x/R -> x/R' and z/R -> (z - z')/R'. Denote 
a ring coordinate counted off the center of a torus cross- 
section (Fig. 1) by x' — R' — R, and therefore, its radius 
is R' = R + a/. We may conveniently introduce the dimen- 
sionless coordinates rj = x'/R, f' = z'/R and p = x/R, 
£ = z/R that will result in an expression for the dimension- 
less potential of the form 



(l + r,')-m r 



■ K(m T ) 



where 



4p- (1 + v') 



(i + v' + p) 2 + (C-C) 2 



(5) 



(6) 



From a condition of the torus homogeneity, mass-to-radius 
ratios for the central and arbitrary component rings are the 
same, and thus, M c = M r /(1 + rj'). Expression for the po- 
tential of the component ring is then 



<Pr(p,C;ri',C') = 



GM C 
■kR 



(7) 



where <f> r is determined by expressions (|6}. Due to ad- 
ditivity, the torus potential can be represented as the in- 
tegral over potentials of the component rings. To do this, 
we replace a discrete mass of the ring M c in ([Til by a 
differential dM, which for the homogeneous torus equals 
dM — -^ijdrj'd^', where M is a total mass of the torus 

7rr Q 

equaling to a sum of masses of the component rings, and 
ro = Ro/R is a dimensionless minor radius of the torus (a 
geometrical parameter) . Then, the potential of the homoge- 
neous circular torus takes the form 



<Ptorus(p, C) 



GM 
n 2 Rrl 



_<Pr(p,C;v',(')dv'd(' 

(8) 

This integral expression for the torus potential is valid for 
both the inner and outer points. The validity of expression 
JH} is confirmed by calculation of the potential made by di- 
rect integration over the torus volume. Hereafter, in analyz- 
ing approximate expressions, we will use the term "exact" 
for the values of potential obtained from the integral for- 
mula (JSj) . In Fig. 2, dependencies of the torus potential on 
the radial coordinate are presented, which were obtained nu- 
merically from formula {SJ for tori with different values of 
the geometrical parameter ro. The potential curves for all 
values of ro are seen to be inscribed into the potential curve 
of the infinitely thin ring of the same mass and radius, lo- 
cated in the torus symmetry plane. The potential curve to 
the right of the torus surface (p > I + ro) virtually coin- 
cides with the potential curve of the ring, while to the left 
(p < 1 + ro), it passes lower and differs by a quantity that 
depends on ro (see section 3). In Fig. 3 the dependencies of 



Here, the coordinates of a point where the ring intersects the 
plane of the torus cross-section are meant as the ring coordinates. 
Position of this point is determined by the ring radius R' and by 
the ring distance from the torus plane of symmetry z' . 
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Figure 2. Dependence of the potential on the radial coordinate 
for C = 0, for the tori with various values of geometrical parame- 
ter: 1) r = 0.05, 2) r = 0.1, 3) r = 0.2, 4) r = 0.3, 5) r = 0.4. 
The potential of an infinitely thin ring {1} with mass M equal to 
the mass of the torus is shown by the dashed line. In this plot 
and in all the subsequent figures, M = 1, R = 1, G = 1. 
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Figure 3. Potential of the torus with ro = 0.3 as a function of 
the radial coordinate for various values of f = 0, 0.1, 0.2, 0.3, 0.4. 



the torus potential on the radial coordinate are presented, 
which were calculated from expression © for different val- 
ues of (. 

Note, that in contrast to the work by Kondratyev (2007, 
p. 196, expression (7.26)), where the torus potential is ex- 
pressed only through a single integration of the elliptical 
integrals of all the three kinds, the torus potential (|8} in 
our work is expressed by double integration of the ellipti- 
cal integral of the first kind. However, further analysis of 
this expression for the torus potential (|8]) allows us to ob- 



tain approximations that are physically understandable and 
enable solving practical astrophysical tasks which need mul- 
tiple calculations of the gravitational potential of the torus. 

For further analysis of the torus potential, we define the 
inner region as the volume bounded by the torus surface (in- 
side the torus body) and outer region as the region outside 
this surface. 



3 TORUS POTENTIAL IN THE OUTER 
REGION 

It is seen from Fig. 2 that the outer potential of the torus 
can be approximately represented by the potential of an in- 
finitely thin ring of the same mass up to torus surface. For 
p — > 0, the values of the torus potential and potential of 
the infinitely thin ring differ by a quantity that depends 
on a geometric parameter ro, that is especially evident for 
a thick torus (ro > 0.5). Find a relationship between the 
outer potential of the torus and the potential of a ring of 
the same mass, that is, derive an approximate expression 
for torus potential in the outer region, where a condition 
(p - l) 2 + < 2 > r\ holds. Within this region, the integrand 
4> r (p, C; rf , C') in (O does not have singularities for all ?/,£', 
therefore, it can be expanded as the Maclaurin series in pow- 
ers of J?',C' hi the vicinity of a point rf = f = 0. Since the 
integrals in symmetrical limits from the series terms that 
contain cross derivatives and derivatives of the odd orders 
are equal to zero, only summands with the even orders re- 
main in the expansion. With the quadratic terms of the se- 
ries being restricted, the potential of the component ring 
is: 



<Mp,C;»/,C') ~ Mp,C)+ 



1 2 <Pr 



2 drf- 



is , 1 
V + — 

v'=o 2 
C'=o 



d 2 <t> r 



d('- 



7J =0 

C'=0 



Substituting <(9j into (|8]), we will have after integration 

d 2 & 



(9) 



<fltor 



GM 
irR ' 



1 + 



2(j> c 



dr) 1 '- 



v '=o d(> 2 
C'=o 



7j'=0 

C'=0. 



(10) 

Ultimately, the approximate expression for the torus poten- 
tial in the outer region (p — l) 2 + £ 2 > r 2 , has a form: 



GM 
nR ' 



(ii) 



where 4> c = J^K(rri) is a dimensionless potential of the 
central ring ((2j), and 



s(p,0 



E(m) 

{p + l) 2 + C 2 ' K(m) 



p 2 + C 2 



(12) 



E(m) = /q 2 d/3\/l — msin 2 /3 is the complete elliptical in- 
tegral of the second kind. We may conveniently proceed to 



a new variable 77 = p — 1 that allows expression (I12|l to be 
represented as 



S(r?,C) = 



r/ 2 + C 2 + 2?7 E(m) 



where 



(r, + 2) 2 + C 2 K(m)' 

! ^ + 1 

(r/ + 2) 2 +C 2 - 



(13) 
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Expression (JTTJ for the torus potential (we will further call it 
the S-approximation) , with (|12l) or (|13[) taken into account, 
represents the torus potential accurately enough in the outer 
region rf + C, 2 > r ( 2 (Fig. 2). Since | S |< 1 the second mul- 
tiplier in is a slowly varying function in p and Let us 
simplify the expression JTT} replacing the second multiplier 
by its asymptotic approximations. 

In the first case, p — > corresponding to r\ — > —1, the 
parameter m — > and E(m)/K(m) — > 1, therefore, S — > 
(£ 2 — 1)/(C 2 + !)• The expression for the torus potential in 
this case is 



<Ptorus(p, (;r ) 



GM 
nR 



MpX) 



1 



16 16 C 2 + 1 



(14) 

Since the dimensionless potential of the ring j2} at the sym- 
metry axis is 4>c = , = , we get for the torus 



<Ptorus(0, C; ?"o) 

and for £ — 0, 



GM 



R yr+c 2 " 



1 r,U 2 -l 
16 16 C 2 + 1 



¥?tor«s(0, 0; r ) « — — I 1 - — 



(15) 
(16) 



The second summand GM/R ■ r /8 in (|16[) describes dis- 
placement of the torus potential at the symmetry axis as 
compared to the potential of an infinitely thin ring (Fig. 2). 

In the second case, at large n, the parameter m — > 0, 
and S — > 1 in (|11[) . therefore: 



GiW 

<y3i O m S (p,C;r ) « — 5"0c(p,C), 
7m 



(17) 



that is, the torus potential is equal to the potential of the 
infinitely thin ring with the same mass M and radius R in 
this case. It is seen from Fig. 4 that the S-approximation for 
the torus outer potential is applicable up to the torus 
surface (upper curves). Indeed, in the region p < 1 — ro, 
difference between the potential obtained from the integral 
expression ^ and its value taken from the S-approximation 
reaches maximum near the torus surface and does not exceed 
0.2% for ro = 0.5. The difference remains small even for 
a thick torus: it does not exceed 1.5% for ro = 0.9. For 
£ = ro all the points are outer, and the curves for the exact 
potential and S-approximation virtually coincide (deviation 
is less than 0.1%). 

Note that asymptotics of the ^-approximation for the 
outer potential (|15p and (|17|l also describe the torus poten- 
tial well enough (dotted line in Fig. 4). Thus, for |£| < ro, 
the approximation (|15|l can be used to estimate the potential 
inside the region bounded by a cylinder with radius p — ro, 
while the approximation (|17[) is applicable outside the re- 
gion bounded by a cylinder with radius p + ro. At |£j ^> 1, 
expression (|15|l tends to HZ}, and expression for potential 
of the infinitely thin ring ((2} can be used within the whole 
outer region to approximately evaluate the torus potential. 

Therefore, the outer potential of the torus can be repre- 
sented with good accuracy by a potential of an infinitely thin 
ring of the same mass. The dependence of the geometrical 
parameter ro appears only in the torus hole; it is taken into 
account in the "shifted" potential of the infinitely thin ring 
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Figure 4. Dependence of the torus potential with ro = 0.5 on p 
for C = (upper curves) and £ = 0.5 (lower curves). Solid lines 
show the potential calculated with the exact formula j5). The S- 
approximation lllll l of the potential is shown in dashed lines, and 
dotted lines represent the limiting cases of the S-approximation: 
the potential curve for the infinitely thin ring H17I I is to the right 
of the torus cross-section, and the curve representing a "shifted" 
potential of the infinitely thin ring (114 6 is to the left. The bound- 
aries of the torus cross-section are dotted with the vertical lines. 



These approximations are valid up to the surface of the 
torus. 



4 TORUS POTENTIAL IN THE INNER 
REGION 

To analyze the inner potential of the torus, it is convenient 
to select the origin of a coordinate system in the center of 
the torus cross-section (Fig. 5) . Then, the dimensionless po- 
tential of the central ring takes a form: 



Mv,0 = 



1+7) 



K(r 



(18) 



where 



4(1 + 7)) 

771 — — 

(2 + r?) 2 + C 2 ' 

Consider the potential of the central ring (|18|) in the vicinity 
of n — > 0, — > 0, that corresponds to m -> 1. In this case, 
the elliptical integral in (jTHJ) can be expanded in terms of 
a small parameter mi = 1 — m. With the series clipped by 
two terms, we will have: 



K(mi) w In- 



4 1 

— = + -mi m — 

/mi 4 e x 



(19) 



3 There is some analogy with the known result: the outer poten- 
tial of a solid sphere of mass M is the same as that generated 
by a point mass M located at the sphere's center. Note, how- 
ever, that torus has another system of equigravitating elements 
(Kondratyev, 2007). 



Gravitational potential of a torus 5 




Figure 5. The scheme of torus cross-section. 



of the component ring: 

<M»7> C; »/> C) ~ f 1 ' f 2 - 



(25) 



where 



In consideration of the inner potential of the torus, rewrite 
expression ([8} in the polar coordinates (Fig. 5): 

GM f r ° f 27r 
<Ptorus{r,9;r ) = / / cj> r (r,0;r ,9 )r dr'dO' ', 

T 2 i»o Jo Jo 

(26) 

where coordinates of the component ring are rf = r'cos#', 
f = r'sinO', and coordinates of a point P are 77 = rcos6>, 
C = r sin 6». Substitute (J23J) and ([24} into and after 

multiplying, restrict ourselves by the terms quadratic in 77, 
£ and 77', f'. Then, after integration of ((26} , we obtain the 
approximate expression for the inner potential of the torus: 

<Ptorus(v,C r o) « 7^ \c + a-irj + a 2 -n 2 + fe 2 < 2 ] , (27) 
where 

fc=^, c= l + 2fc 2 -21nfc + 8fc 2 lnfc, 5i = l + lnfc, 



mi 



T- 2 + 4(1 + 77) ' 



r 2 = 77 2 + f 2 . The approximate formula for the ring potential 
expressed through the parameter mi is then: 



Mv,0 



In- 



/mi 



1 1 4 
—mi in — ■ 

4 eJm\ 



(20) 



Passage to the potential of an arbitrary component ring is 
fulfilled by substitutions 1 + 77 — > (1 + rj)/(l + 77') and £ — ► 
(f — C')/(l + ^Oi which results in an expression 



2(1 + 7/) 



In 



1 + 



'"1 



mi 
~4~ 



(21) 

where m[ = (r - r') 2 /g, g = (r - r') 2 + 4(1 + 77) (1 + 7/). A 
summand (r — r') 2 = (77 — ?7') 2 + (£ — £') 2 is a square of the 
distance between the component ring and a point P (Fig. 5). 
Expansion (|21l) is valid for m'x — > 0, therefore, (r — r') 2 <C 1. 

Confine ourselves by the case of a thin torus (ro 
Then (r-r') 2 < 4(1 + rj) (1 + 77') and ( 1 + 77) ( 1 + r?') w 1, and 
the first multiplier in ()2ip can be written to the second-order 
terms as: 



h 



2(l + r?Q 



I + 77' 
l + r? 



(22) 



After expanding the square root in (|22|i in powers of 77 and 
77', we obtain 



h 



77' 77 7777' 77' 3 2 



1 If '\2 

l-g(r-r) 



(23) 

Similarly, the second multiplier (in square brackets) in ex- 
pression (|21[) can be written to the terms quadratic in coor- 
dinates as: 



h 



In 



1 + 



'"1 



'"1 

T 



1 



-fin- 



8p 



(24) 



r — r' I 16 I r — r' 

Thus, we obtain the following expression for the potential 



1 



a 2 



(8k) 2 



4fc 2 (ll + lOlnfe), 



(8fc) 



' +4fc 2 (3 + 21nfc). 



b, = 



The first summand in (|27[) is the value of the torus potential 
in the center of the torus cross-section: c = <j>torus{0, 0; ro). 
To further analyze the inner potential, it is convenient to 
transfer to a coordinate system normalized to the geomet- 
rical parameter of the torus ro- Then the series coefficients 
will transform to the form: 

ai = 8Jfe(l + lnfc), a2 = -1 - 4fc 2 (ll + lOlnfc), 

b 2 = -l + 4fc 2 (3 + 21nfc). 

The expression for the torus potential (|27[1 can be written 
as 



<fitorus(r],C\ r o) 



GM 
2-kR 



ai h a- 2 

ro 



ro 



ro 



(28) 

It follows from (|28|) that the maximal value of the poten- 
tial reaches at a point rj m ax ~ — (airo)/(2a2), C — 0| while 
equipotential lines are ellipses with their centers displaced an 
amount rjmax with respect to the center of the torus cross- 
section, and a ratio of semiaxes of the ellipses is \fb2~Jo2- 
Note, that location of the potential maximum 77 = rjmax, 
£ = corresponds to the weightlessness point, where the re- 
sultant of all the forces affecting a particle inside the torus 
equals zero. In such an approximation, components of the 
force inside the torus depend on the coordinates linearly at 
that. In Fig. 6, the curves of the inner potential in the coor- 
dinate system normalized to ro are presented for three values 
of r . 

Though we confined ourselves to the case of a thin torus, 
the curves of potential taken from expression (|28[) are well 
consistent with the curves for the exact potential © up to 
ro = 0.5, where the deviation is maximal near the torus 
surface and is of the order of 2% (Fig. 6). The value of 
potential in the center of the torus cross-section (a constant 
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Figure 6. Dependence of the inner potential on the normalized 
coordinate n/ro at f = for various values of the geometric pa- 
rameter: ro = 0.1,0.2,0.5. Solid curves represent dependencies 
of the potential on the distance to the center of the torus cross- 
section, which were calculated with the exact formula {§}. Dashed 
lines are the potential curves taken with the approximate expres- 
sion for the inner potential of the torus (1 2S D - 



c = 4>torus(0, 0; ro) in (1271) ). also coincides with its exact 
value. 

It is of interest to investigate the solutions obtained for 
limiting cases. Indeed, the case ro = Ro/R corresponds 
to two limiting passages to an infinitely thin ring (R is fixed 
while Ro — > 0) and to a cylinder, when Ro is fixed and 
R — > oo. Dwell on the limiting passage to the cylinder po- 
tential. At ro — ¥ 0, the coefficients c — > 1 — 2 In k, a\ — > 0, 
and coefficients a2,&2 — > — 1, and expression (1281) takes the 
following form in this case: 



<ptaru S {r],C;ro -> 0) w — - 



i_ 21n !£_(£. 

ro 



(29) 



where r = r) + f . It is known that the inner potential 
of a circular cylinder with the length 2H, much larger than 
the radius Ro of its cross-section (Kondratyev 2007) , has the 
form: 



GM 
~2H 



2 In 



2Veff 
Ro 



ro 



(30) 



After a formal substitution 2H — 2-kR in (|30[1 , (the cylinder 
length is equal to the length of the central ring), we get an 
expression: 



GM 
2nR 



\2ty 



(31) 



Expression (|3ip coincides with Q29p to a constant Q- 

The quadratic dependence on r for the inner potential 

4 Difference in the constant may be caused by the curvature of 
the torus surface. 



of a thin torus can be also derived in the case, when the 
minor radius Ro — > 0. The outer potential of the torus was 
shown in section 3 to be approximately equal to the potential 
of an infinitely thin ring of the same mass and radius. In 
this case, the smaller is the torus geometrical parameter ro, 
the more accurate is this approximation. Therefore, at rj 2 + 
C 2 > r o — > 0, the outer potential of the torus tends to the 
potential of an infinitely thin ring. In this case, rj, £ — > 0, and 
thus, the elliptical integral in the expression for an infinitely 
thin ring Q can be expanded in the vicinity of m — > 1. If 
we confine ourselves to the first term of the expansion, we 
get an approximate expression for the potential of a central 
infinitely thin ring 

GM 



2-kR 



(-ln(j7 2 + C' ! ) + 21nf 



(32) 



which remains valid for the outer potential of the thin torus 
as well. It should be noted that there is no dependence on ro, 
because in such an approximation, all thin tori with the same 
masses and major radii are equigravitating for the outer 
potential. The derivatives of the outer potential of the thin 
torus in rj, £ are then: 

d^ r ._GM_v_ dip c ^ GM C 
dr] ttR r 2 ' irR r 2 

and take the following forms at the torus surface (rj 2 + ( 2 — 
ro 2 ): 



d(p c 
drj 



GM 
nRro 



difc 



GM 
nRro 



It is the linear dependence of the force on coordinates rj, £ 
that satisfies such boundary conditions. Thus, the inner po- 
tential of the thin torus can be represented to the integration 
constant in the form: 



<Ptorus(r], C; r 0) 



GM 
2-kR 



n 2 + C 2 



+ c(r ] 



(33) 



Equating ([32} with (|33[) at the torus surface, we obtain the 
expression for the constant c(ro) = — 21n(ro/8) + 1 that 
coincides with expression (|27[) obtained above at ro <C 1. 

It becomes evident from analysis of the inner potential 
for the two limiting cases (Ro — > and R — > oo) that the 
first summand in coefficients 02, 62 of the power series (|28[l 
represents properties of the inner potential of a cylinder. 
With the cylinder potential separated, the inner potential 
of the torus (|28p can be written as: 



<Ptorus(r], = <ficyl(r) + Ificurv^, C)> 



where 



^Pcurv 



GM 
2-kR 



- ai 



+ C a 



+ c b 



(34) 



(35) 



21n— + 2/c (l + 41nfc), 

27T 



1 



«2, 



I+62 



The second summand <p C urv(rj, C)i which we will call a po- 
tential of curvature, implies curvature of the torus surface. 
Indeed, all the coefficients of the series (|35jl tend to zero in 
the limiting passage to the cylinder (ro — > 0), and <p C urv — > 0. 
Therefore, the inner potential of the torus can be represented 
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as a sum of the cylinder potential and a term comprising a 
geometrical curvature of the torus surface. 



5 SEWING TOGETHER THE INNER AND 
OUTER POTENTIALS AT THE TORUS 
SURFACE 

In the previous sections, we derived approximate expres- 
sions for the torus potential in the outer (n 2 + £ 2 > Tq) 
and inner (r/ 2 + £ 2 < r 2 ,) regions. It has been shown also 
that the inner potential of the torus can be represented by 
a series in powers of n/ro and C/ r o, and the constant, linear 
and quadratic terms of the series were determined analyti- 
cally. To find a larger number of the series terms sufficient 
to represent the inner potential accurately enough, and to 
obtain a continuous approximate solution for the potential 
and its derivatives in the whole region that would satisfy 
the boundary conditions at the surface, we will act in the 
following way. Represent the inner potential of the torus as 
a power series [ff 

M.C;n>) = ^fc(ro)+^ai(r )^) + 

+ EEm™>(^)'+Em„)(^), 

(36) 

where c(ro), aj(ro), bj(ro), tij(ro) are unknown coefficients. 
Note, that the series (|36|) contains only the terms with even 
powers of £, because the torus potential is symmetric in £. 

Suppose that we have an analytical expression for the 
torus potential ^(n, f; ro) at its surface (n 2 + £ 2 = r§). Also, 
write down the inner potential of the torus (|36[) at its surface 
(n = ro cos 9 and £ = ro sin 9) : 

<j)(0, r ) = i ( c + ^2 a.i cos 1 9 + 

J2bjsm'e). (37) 

i=l 3 = 1 3=1 / 

From conditions of equality of the inner and outer potential 
and its derivatives in coordinates at the torus surface for 
several angles 9k , we obtain a system of 3fc linear equations 
to determine coefficients c, a%, bj, tiji 

' c+ E o< cos 4 <9 fe + E E *ij cos 4 6> fc sin J ' 9 k + 

i — l i — l j=l 

+ E bjSin'Ok =2TT^(0 k ,r o ) 

3=1 

E i • <N cos 1-1 Sfc + E E * ' *»i cos 1-1 # fe sin J # fc = 

\ i = l i = l 3=1 

= 2^ro^*(e fejro ) 

E j • >' sin^ 1 6> fc + E E 3 ■ Uj cos 1 9, sin 3 '" 1 fl fc = 

3=1 i=l j = l 

, =27rro^*(0fc,r o ) 

(38) 

5 We consider a dimensionless potential here. To pass to the di- 
mensional case, J 36 1) must be multiplied by GM/R. 
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Figure 7. Dependence for the torus potential on p for ro = 0.5 
(£ = 0) in the whole region: the potential curve calculated from 
the exact expression (JSJ is shown by a solid line; dashed line 
demonstrates a result of sewing together the S-approximation in 
the outer region with the inner potential represented by the 
power scries II37I I with the coefficients found from the sewing con- 
dition l |39) . 

Thus, if we had the analytic solution for the outer poten- 
tial of the torus, we could obtain an exact expression for 
the inner potential as an infinite series in powers of cos 9k, 
sin 9k, using the boundary conditions and solving the sys- 
tem of equations (|38[) . Since there is no analytic expression 
for the outer potential, we can use the above approximate 
expression (If I ll for the torus potential in the outer region 
(the S-approximation), and introduce designations: 

* = E [<t>in(9k, ro) — 4>out{9k, ro)] 2 
fc 

*i =E [fc»(0k,ro) - 4w(0fc,r o )]) 2 

$2 = E [<Pin(e k ,r ) - 0out(8 k ,r o )]y , 

where <f>i n , cf> m t are solutions for the inner (|37(l and outer 
HI 1 1 potentials at the torus boundary, respectively. The un- 
known coefficients of the series can be then determined from 
a condition of the minimal value of a functional: 

F = $ + $i + $2 -> min. (39) 

The functional (|39[) was minimized with the least squares 
method, and coefficients of the series (|37(l were determined 
up to the 4-th power. The coefficients of the series are pre- 
sented in Appendix (Table Al). In Fig. 7 dependence of 
the potential on the radial coordinate from the exact ex- 
pression ([8]) is presented for the entire region, as well as 
its approximate solution obtained by sewing together the 
S-approximation and the inner potential ()37[) . Though 
the approximate solutions were obtained assuming that the 
torus is thin ro <C 1, we see that the exact (JHJ) and ap- 
proximate solutions are consistent even for the torus with 
ro — 0.5. In Fig. 8, the equipotential curves on the plane of 
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Figure 8. Equipotcntial curves for the torus with rrj = 0.5. Solid 
lines are those calculated from the exact expression for the poten- 
tial; dashed lines are those for the approximate formulas for the 
outer and inner potentials (see the figure 7 caption). The torus 
body cross-section is indicated as a gray circle. 

the torus cross-section are shown, where a good agreement 
for all values of p, £ is seen as well. 



6 CONCLUSIONS 

In the present work, the gravitational potential of a homo- 
geneous circular torus is investigated in details. An inte- 
gral expression for its potential that is valid for an arbitrary 
point is obtained by composing the torus of infinitely thin 
rings. This approach has made it possible to find an ap- 
proximate expression for the outer potential of the torus 
(S-approximation), that has a sufficiently simple form. It is 
shown that the outer potential of the torus can be repre- 
sented with good accuracy by a potential of an infinitely 
thin ring of the same mass. The dependence of the geomet- 
rical parameter ro appears only in the torus hole; it is taken 
into account in the "shifted" potential of the infinitely thin 
ring. These approximations are valid up to the surface of the 
torus. 

For the inner potential, an approximate expression is 
found in the form of a power series to the second-order terms, 
where the coefficients depend only on the geometric parame- 
ter ro. Expressions for the potential in the center of the torus 
cross-section and for coordinates of the potential maximum 
are obtained, and the limiting passage to a cylinder poten- 
tial is considered. It is shown that the inner potential of the 
torus can be represented as a sum of the cylinder potential 
and a term comprising a geometrical curvature of the torus 
surface. A method for determining the torus potential over 
the whole region is proposed that implies sewing together 
at the surface of the outer potential (S-approximation) with 
the inner potential represented by the power series. This 
method provided a continuous approximate solution for the 
potential and its derivatives, working throughout the region. 

Surely, matter distribution within a torus is inhomo- 
geneous in actual astrophysical objects, and a torus cross- 
section may differ from a circular one. Therefore, it is further 
interesting to account for inhomogeneity of matter distribu- 
tion inside a torus, for difference of the torus cross-section 
from a circular form, and so on. 
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APPENDIX A: COEFFICIENTS OF THE 
POWER SERIES FOR INNER POTENTIAL OF 
A TORUS 

In Table Al, coefficients of the power series (up to the 4-th 
power) for the inner potential of the torus, are presented, 
which were calculated from the sewing condition for the 
torus with various values of the geometrical parameter ro. 
The analytic expression (|28[) was used to determine the zero- 
th coefficient c of the series. In Fig. Al, the linear (ai) and 
quadratic (d2, 62) coefficients of the power series as functions 
of ro obtained analytically from (|27[) are shown by solid lines; 
the dots are the proper values of these coefficients obtained 
from the condition of sewing (see Table Al). The values of 
the analytic coefficients are seen to coincide up to ro = 0.5 
with their values obtained independently with the method 
of sewing from (1391) . 
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Table Al. Coefficients of the power series for the inner potential of the torus for various values of rrj, 
obtained with the method of sewing. 



Coeff. Geometrical parameter rrj 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 




0.7 


0.8 


0.9 


Ol 


-0.33798 


-0.53651 


-0.68154 


-0.79129 


-0.87439 


-0.93587 


-0 


.97906 


-1.00628 


-1.01928 


0.2 


-0.98002 


-0.93543 


-0.87773 


-0.81171 


-0.74107 


-0.66865 


-0 


.59677 


-0.52739 


-0.46224 


f>2 


-1.00411 


-1.01086 


-1.01970 


-1.02892 


-1.03759 


-1.04495 


-1 


.05030 


-1.05299 


-1.05237 


a-A 


0.02392 


0.04364 


0.05781 


0.06608 


0.06853 


0.06550 





.05753 


0.04525 


0.02938 


*12 


0.02550 


0.05329 


0.08404 


0.11791 


0.15454 


0.19323 





.23295 


0.27238 


0.30991 


a 4 


-0.00182 


-0.00785 


-0.01610 


-0.02576 


-0.03580 


-0.04535 


-0 


.05371 


-0.06036 


-0.06495 


f>4 


0.00061 


0.00131 


0.00308 


0.00570 


0.00922 


0.01362 





.01880 


0.02453 


0.03045 


*22 


-0.00122 


-0.00812 


-0.01948 


-0.03681 


-0.06076 


-0.09157 


-0 


.12900 


-0.17213 


-0.21931 



